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ABSTRACT 

We have observed six high-mass star-forming regions in the ^113/2, J ~ 7/2 Hnes of OH using 
the GBT in order to investigate whether the magnetic field, and hence the density, measured in 
absorption differs from that implied by maser Zeeman splitting. We detect absorption in both 
the f344f and 13434 MHz main lines in all six sources. Zeeman splitting in the = 3+ 3^ 
absorption line in W3(0H) implies a line-of-sight magnetic field strength of 3.0 ± 0.3 mG. This 
is significantly less than full magnetic field strengths detected from OH maser Zeeman splitting, 
suggesting that OH maser regions may be denser than the non-masing OH material by a factor 
of several. Zeeman splitting is not detected in other sources, but we are able to place upper 
limits on _B|| of 1.2 mG in GIO. 624— 0.385 and 2.9 mG in K3— 50. These results are consistent 
with a density enhancement of the masers, but other explanations for the lower magnetic field in 
absorption compared to maser emission are possible for these two sources. Absorption in one or 
both of the 13442 and 13433 MHz satellite lines is also seen in four sources. This is the very first 
detection of the ^113/2, J = 7/2 satellite lines. Ratios of satellite-line to main-line absorption 
suggest enhancement of the satellite lines from local thermodynamic equilibrium values. Masers 
are seen in the F = 4+ ^ 4~ and 3+ — > 3~ transitions of W3(0H) and the 4+ 4~ transition 
of ON 1. A previously undetected 4+ 4^ maser is seen near —44.85 km s^^ in W3(0H). 

Subject headings: masers — ISM: magnetic fields — H II regions — radio lines: ISM — stars: 
formation — ISM: molecules 
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1. Introduction 

Hydroxyl (OH) masers are found in higli-mass star-forming regions. They trace the local velocity and 
magnetic field and therefore provide clues to understanding the physical conditions of material surrounding 
newly-formed high-mass stars. But an unanswered question is whether the physical conditions in masing 
regions are representative of the surrounding material. In order for OH masing to occur, a large column 
density of OH must have velocity coherence such that the velocity gradient along the amplification path 
of the masing clump does not exceed the masing linewidth. In principle, sufficient column density can be 
achieved in two ways. If a masing clump is much denser than the ambient gas, the total OH column density 
can be large even if the physical extent of the clump is not. Alternatively, masing may occur along favored 
paths of velocity coherence even in a medium of homogeneous density. 

While it is uncertain which of these two scenarios is prevalent, two pieces of evidence suggest that 
density enhancements may not be necessary for OH masing to occur. First, 13434 MHz OH absorption 
in W3(0H) shows Zeeman splitting indicating a line-of-sight magnetic field strength of 3.1 mG (Glisten, 
Fiebig, & Uchida 1994). VLBI measurements of 1665 and 1667 MHz OH masers in W3(0H) imply similar 
magnetic field strengths (e.g., Bloemhof, Reid, & Moran 1992). The magnetic field strength of collapsing (or 
collapsed) material scales as the density n'^, where n w 0.5, a result supported both by theoretical modelling 
(Mouschovias 1976; Fiedler & Mouschovias 1993) and observations of molecular clouds in various stages 
of collapse (e.g., Crutcher 1991). Since the OH absorption and maser magnetic fields are comparable, this 
suggests that the density at masing sites is similar to that of the ambient cloud of OH. Second, ammonia 
(NH3) observations of W3(0H) indicate that the density of material in the clumps of maser emission is 
roughly the same (to within a factor of two) as the density of the interclump material (Reid, Myers, & 
Bieging 1987). Since the velocity and extent of NH3 absorption is similar to that of the OH emission in 
W3(0H), it is a reasonable assumption that NH3 and OH exist in the same cloud of material. Finally, 
Cesaroni & Walmsley (1991), when modeling multi-transition OH observations of W3(0H), find that maser 
emission in certain lines and absorption in others can be explained for the same range of densities, between 
10^ and a few times lO'' cm~^ (for a temperature of 150 K). 

Is W3(0H) a special case, or is it representative of all interstellar OH maser sources? In particular, is 
any density enhancement required for the onset of OH masing? These are questions that motivated us to 
observe a wider range of sources than Glisten, Fiebig, & Uchida (1994) with the high sensitivity that the 
GBT can afford. 



2. Observations 

The observations were performed on 2004 April 11 and 12 using the National Radio Astronomy Ob- 
servatory's'^ Robert C. Byrd Green Bank Telescope (GBT) in Green Bank, WV. The GBT has an effective 
diameter of 100 m. Observations were taken in both circular polarizations with the Gregorian focus Ku- 
band receiver. The GBT Spectrometer was configured in 9-level mode to provide 8192 uniform-weighted 
spectral channels in eac;h of two IFs covering a bandwidth of 12.5 MHz each, centered on 13441.4173 and 
13434.6374 MHz, Doppler-shifted to the LSR velocity of each source. These correspond to the frequencies 
of the F = 4+ — > 4~ and F = 3+ — > 3~ main-line ^H3/2, J = 7/2 transitions, respectively (Destombes et al. 



^Thc National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative 
agreement by Associated Universities, Inc. 
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1977). Each spectral channel had a resolution of 1.5 kHz, corresponding to a resolution of 0.034 km s ^ in 
velocity space. Thus, the 13433.982 MHz (i^ = 4+ ^ 3") and 13442.072 MHz (i^^ = 3+ ^ 4") satellite lines 
were also in the observed frequency range. At these frequencies, the FWHM beamwidth of the GBT is 55" 
and the gain is 1.7 K/Jy. 

Double beam switching was employed such that the source appeared in each of the dual beams alternately 
for a period of two minutes. The beams were separated by 330" in azimuth. Dynamic pointing and focusing 
were used, and the pointing and focus were checked hourly as well as each time the telescope was pointed 
at a source. 

Variable rain, heavy at times, fell throughout the data collection period. During the brief interludes 
without precipitation, the (non-zenith) system temperature approached being receiver-limited at about 30 K. 
At times, the system temperature was over 100 K, due to heavy precipitation and low elevation angles. 

Our sources were chosen according to two criteria. First, they must contain a strong background H II 
region so that absorption might be observed. Second, since it was our intent to compare the magnetic field 
strength in the OH gas seen in absorption with that seen in maser emission, we sought sources in which the 
masers indicated a uniform field distribution and eliminated any with a reversal of the line-of-sight direction 
of the magnetic field as determined from ground-state OH maser Zeeman splitting (Fish 2004). Only about 
10 sources observable at the latitude of the GBT meet these criteria. We observed the six sources listed in 
Table 1. 



The masers we detected are shown in Figures 1 to 4 and discussed in §3.1. The absorption features we 
detected arc shown in Figures 5 to 21 and discussed in §3.2. The absorption spectra have been Hanning 
weighted for clarity, although the analyses are based on the uniform-weighted data. The antenna tempera- 
tures are related to the flux density hy S = 2kTA/A, where A is the effective collecting area of the telescope. 
We have used the convention that Stokes I = 0.5 (Ta,lcp + Ta,rcp) and Stokes V = 0.5 {Ta,lcp — 2a,rcp)-^ 



2 We define Stokes V as LCP - RCP to be consistent with Gusten, Fiebig, & Uchida (1994). 



3. 



Results 



Table 1. Observed Sources 



Source 



RA(J2000) Dec(J2000) 



Obs. Time 
(min) 




W3(0H) 

GlO.624-0.385 

G28.199-0.048 

W49 

K3-50 

ON 1 



02 27 03.70 
18 10 28.61 

18 42 58.04 

19 10 11.04 

20 01 45.73 
20 10 09.05 



+61 52 25.4 
-19 55 49.7 
-04 13 58.0 
-h09 05 20.2 
+33 32 45.3 
+31 31 35.2 



276 
270 
10 
58 
178 
20 



0.006 
0.008 
0.026 
0.028 
0.010 
0.015 



'^Single-channel RMS noise in Stokes I and V. 
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Thus, Ta = 1 K in both RCP and LCP would result in S' = 0.59 Jy in Stokes I. Zeeman measurements in 
the OH absorption are discussed in §3.3. Remarks about absorption line ratios are presented in §3.4. 



3.1. Maser Emission 

Masers were found in two sources: W3(0H) and ON 1. Parameters of the maser lines are listed in 

Table 2. Magnetic field estimates assume a Zeeman splitting coefficient of 0.f78 km s~^ mG~^ for the 
4+ — > 4~transition and 0.230 km s~^ mG^^ for the 3"*" 3~transition (Giisten, Ficbig, & Uchida 1994). 
Zeeman splitting of OH masers is sensitive to the strength of the full, three-dimensional magnetic field, 
independent of its inclination to the line of sight, when the Zeeman splitting exceeds the linewidth. 

In W3(0H), we find maser lines in both the '^11^/2, J = 7/2, J" = 4+ ^ 4~ and 3+ 3~ transitions. 
We detect seven or eight maser line components in each circular polarization in the F = 4+ ^ 4^ transition, 
compared to the three previously detected with other single dish antennas (Baudry ct al. 1981; Giisten, 
Fiebig, & Uchida 1994; Baudry & Desmurs 2002). Nearly all of these lines can be grouped into Zeeman 
pairs with implied full magnetic field strength ranging from 6.9 to 11.3 mG. This range agrees with previ- 
ous observations by Baudry & Diamond (1998) with the VLBA, whose angular resolution is sufficient to 
unambiguously pair most Zeeman components. The strongest masers we detect have counterparts in Table 
1 of Baudry & Diamond. It is not possible to identify unambiguous counterparts to our weaker masers, 
which may be blends of features resolved at the submilliarcsecond resolution of the VLBA. The detection of 
the LCP and RCP masers centered at —44.85 km is new. There is some ambiguity as to whether the 
LCP maser at —43.66 km s~^ should be paired with the RCP maser at —43.54 or —43.51 km s~^, but the 
former is more likely because the linewidth of the latter is much greater than for the LCP maser. In the 
F = —> 3~ transition, we find one Zeeman pair implying a linc-of-sight magnetic field component of 10.3 
mG, consistent with the detection by Giisten, Fiebig, & Uchida (1994). 

ON 1 is the only other source in which we find "^113/2, J = 7/2 masers. The F = 4+ ^ 4~ maser at 
14 km s~^ was also detected by Baudry & Desmurs (2002), although their velocity resolution was insufficient 
to find measurable Zeeman splitting. We do not see the strong (0.5 Jy) maser at —0.13 km s~^ that they do; 
however, we find evidence of a weak maser at 0.38 km s~^ (LCP) and 0.24 km s^^ (RCP). While these lines 
are weak, we believe they are real. The peaks of the LCP and RCP lines are only 3.8 and 2.8 times the single- 
channel RMS noise of 0.02 K, respectively, but the line widths are 8 times a single channel width. These 
masers can be interpreted either as F = 4+ — > 4~ lines at 0.3 km s~^ or as F = 3+ — *■ 4~ lines at 14.9 km s~^. 
The 1665 MHz (^H3/2! J = 3/2) masers in ON 1 are grouped in two velocity ranges: 1 to 4 km s~-^ and 10 
to 16 km s^^ (Argon, Reid, & Menten 2000); likewise, the 6035 MHz (^H3/2, J = 5/2) masers are seen from 
-1 to 2 km s^i and 13 to 16 km s^^ (Baudry et al. 1997). Thus, either the F = 4+ ^ 4" or F = 3+ ^ 4" 
interpretation of these maser lines is consistent with the ground-state OH maser velocities. 

Nevertheless, two arguments suggest that the weak maser feature in ON 1 is a maser in the F = 4+ — » 4~ 
transition. First, satellite-line emission in interstellar OH masers is generally weak compared to main-line 
masers. This is true in the ^11^/2, J = 3/2 lines, where the 1665 and 1667 MHz masers are typically much 
stronger than 1612 and 1720 MHz masers. In the ^H3/2, J = 5/2 lines, strong masers are found in the 6035 
and 6030 MHz main lines, but 6049 MHz satellite-line emission is weak and rare (Baudry et al. 1997), while 
the 6016 MHz satellite line has not been seen in emission (Baudry et al. 1997; Gardner & Martfn-Pintado 
1983). Indeed, if photon trapping is important, inversion in the 6016 MHz satellite line may be impossible 
(Elitzur 1977). As for the ^H3/2, J = 7/2 lines, while about a dozen masers have been found in the 13441 
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MHz main line transition, the 13434 MHz maser in W3(0H) remains the single known maser in the other 
main line transition, despite two searches encompassing 77 distinct interstellar masing sources (Baudry & 
Desmurs 2002; Caswell 2004). Empirically, this suggests that satellite-line ^113/2, J = 7/2 masers should be 
extremely rare if they exist at all. Second, due to a smaller Zeeman coefficient, the magnetic field strength 
implied by a constant-velocity Zccman splitting is higher for satellite lines than for main lines. If the weak 
maser in ON 1 is actually in the F = 3+ ^ 4~ transition, the implied magnetic field is —14.5 mG. This field 
strength would be much larger than that derived from any other Zeeman pair in ON 1 in the ^113/2 ladder 
(see also Baudry et al. 1997; Fish 2004). 

3.2. Absorption 

Unlike maser emission, absorption is clearly seen in all sources in our sample in both the f = 3+ ^ 3~ 
and F = 4+ ^ 4~ transitions. In some sources, absorption is seen in the F = 3+ — > 4~ and/or F = 4+ — > 3~ 
transitions as well. We believe that this is the first clear detection of absorption in the ^n3/2, J — 7/2 satellite 
lines. Parameters of the absorption lines are given in Table 3. We summarize results for individual sources 
below. 

W3(0H): Absorption is seen in all four hyperfine transitions in W3(0H). The presence of multiple 
strong maser lines (see Figures 1 and 5) in the F = 4+ — > 4~ transition is problematic for fitting Gaussian 
components to the absorption. The data in the velocity range from —45.52 to —41.00 km s~^ have been 
excluded from the fitting of the two Gaussian components listed in Table 3. 

GIO. 624— 0.385: Absorption is seen in the F = 3+ — > 4~ transition as well as the main lines. 

G28. 199—0.048: Absorption is seen in the main lines only. The noise is high due to the short time spent 
on source. 

W49: Very broad absorption is seen in the main lines. This absorption is broader than the frequency 
difference between adjacent main- and satellite-line transitions. It is possible that the absorption marked at 
5.59 km s~^ in the F = 3+ ^ 4^ transition is actually at —9.01 km s~^ in the F = 4+ ^ 4~ transition. 
However, the velocity assuming that this component is due to the F = 3+ ^ 4~ transition is more consistent 
both with the other two transitions and with the H92a/He92a velocity of 6.2 ± 2.3km s~^ for the associated 
H II region (Dc Prcc, Mchringcr, & Goss 1997). 

K3— 50: Absorption is seen in all four hyperfine transitions. 

ON 1: Absorption is seen in the main lines only. While the observations occurred during a period devoid 
of precipitation, the noise is still high due to the short time spent on source. 

3.3. Zeeman Measurements in Absorption 

The magnetic field can be measured by Zeeman splitting in absorption, as well as in maser emission. 
While the full three-dimensional magnetic field strength is obtained from Zeeman splitting of velocity- 
separated maser components, absorption lines are much broader than the velocity separation between LGP 
and RCP components, so only the line-of-sight component of the magnetic field can be measured (e.g., Sault 
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et al. 1990). The Stokes V curve is related to the derivative of the Stokes I curve by the following equation: 

^vM-f^^B,, (1) 

where Tv(i^) and Ti(i^) arc the brightness temperatures of the Stokes V and I spectra, and C = —1.06 x 
10*^ Hz for the = 3+ -> 3" transition and 7.95 x 10^ Hz G"^ for the 4+ 4" transition (Glisten, 
Fiebig, & Uchida 1994). Hence, if the Stokes V spectrum of an absorption component is nonzero, the 
line-of-sight magnetic field can be measured. 

W3(0H) is the only source in which the Stokes V absorption spectrum shows clear evidence of Zeeman 
spHtting. Figure 22 shows the Stokes V spectrum of the F = 3+ ^ 3^ transition of W3(0H). The feature 
near —42.4 km s^^ is due to the maser. In the direction of lower velocity from this feature is an "S-curve." 
Superposed atop the data are two curves corresponding to the derivative of the Gaussian fits to the Stokes I 
curve. The curve labelled "Main" indicates the scaled derivative of the Gaussian centered at —45.03 km s~^ 
(Ustcd in Table 3) , while the curve labelled "Both" indicates the scaled derivative of both this Gaussian and 
the one centered at —47.22 km s~^. A positive magnetic field shifts RCP to higher velocity than LCP; in 
absorption, this corresponds to the positive bump of the Stokes V curve being at higher velocity than the 
negative bump. Note that in emission, a positive magnetic field results in the positive bump of the Stokes V 
curve being at lower velocity than the negative bump, as is seen for the 42.4 km s^^ maser Zeeman pair in 
Figure 22. Thus the absorption and emission magnetic fields in W3(0H) are consistent in sign even though 
the S-curves appear inverted with respect to each other. 

The magnetic field value that provides the best fit to the data is 2.9 ± 0.3 mG for the main feature and 
3.0±0.3 mG for both features combined. This is consistent with the values of 3.1±0.4 mG and 3.2±0.6 mG 
obtained by (Giisten, Fiebig, & Uchida 1994). 

For the other sources, no S-curve is observed in the Stokes V spectra, but we can place upper limits on 
the possible magnetic field strength. From Equation 1, the measurement error of the line-of-sight magnetic 
field component in a single channel is 

In principle, each spectral channel constitutes an independent measurement of ctbh • The overall error of the 
parallel magnetic field strength can be estimated as 

O-Sy .overall = (Xl'^Sll) ' ("^^ 

where the summation is taken over all spectral channels. 

Table 4 gives the 3cr upper limits on magnetic field strengths from absorption measurements. For com- 
parison, the average three-dimensional magnetic field strengths obtained from OH maser Zeeman splitting 
are also provided. Masers from each of the J = 3/2, 5/2, and 7/2 sets of transitions are considered sepa- 
rately. When Zeeman splitting is detected in more than one set of transitions, the minimum and maximum 
average magnetic strengths obtained from OH maser Zeeman splitting in sets of transitions corresponding to 
different values of J arc quoted. The actual range of magnetic field strengths seen in maser Zeeman splitting 
may be greater that that listed in Table 4, but averages seem appropriate for comparison with our absorption 
data, since the GET beamwidth is much greater than the angular extent of any of our sources. With the 
exception of the Zeeman pair centered at —44.85 km s~^, Zeeman pairs in the F = 4+ ^ 4~ transition 
detected in this work are excluded from the calculation of the average magnetic field strength in W3(0H) 
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obtained from J = 7/2 masers, in deference to the much higher angular resolution obtained by Baudry & 
Diamond (1998). Note that the magnetic fields derived from the Zeeman pairs we detect agree with those 
obtained by Baudry & Diamond in sign and magnitude to better than 1 mG. 

In addition to obtaining a positive field measurement for W3(0H), we are able to place upper limits of 
several milligauss on B|| for two other sources: K3— 50 and GIO. 624— 0.385. In K3— 50 the magntitude of 
the line-of-sight component deduced from the Zeeman splitting of the main absorption feature is less than 
2.9 mG. Measurements of OH masers give three-dimensional magnetic field strengths of —2.6 to —7.5 mG 
in the ^Bg/s, J = 3/2 hues (Fish 2004) and -5.3 to -9.1 mG in the 6035 MHz J = 5/2 line (Baudry et al. 
1997). 

Our 3a limit of 1.2 mG for GIO. 624— 0.385 improves by a factor of 2 on previous observations in the OH 

^H3/2, J = 7/2 F = 3+ ^ 3^ line and is comparable to results obtained from observations in SO absorption 
(Uchida, Fiebig, & Giistcn 2002). In GlO.624-0.385, VLBA observations of the 1667 MHz J = 3/2 line 
show one Zeeman pair with a magnetic field of —6.0 mG (Fish 2004). 

Four possibilities, alone or in combination, can explain the non-detection of a Zeeman pattern in the 
absorption in K3— 50 and GIO. 624— 0.385. (1) A reversal of the line-of-sight direction of the magnetic field 
exists across each source. (2) The magnetic field in K3— 50 and GIO. 624— 0.385 is inclined at a large angle 
relative to the line of sight. (3) The average magnetic field strength of the material sampled by ^11^^/2, J — 7/2 
OH is smaller than that sampled byJ = 5/2orJ = 3/2 OH. (4) Masers are on average denser than the 
surrounding OH material, and thus the magnetic field strength is greater at masing sites than for the OH 
as a whole. 

The first explanation is possible for both these sources. While five "^^3/2, J = 3/2 maser Zeeman pairs 
were found in VLBA imaging of K3— 50, they are all located on the periphery of the H II region to the north 
and east. It is possible that the line-of-sight direction of the magnetic field reverses on the south or west 
side of the source. Even if there is no reversal in K3— 50, four of the five Zeeman pairs found in the VLBA 
images of the 1665 and 1667 MHz maser emission imply full (unprojected) magnetic field strengths less than 
the dctectability limits in Table 4 (Fish 2004). 

Only one ground-state OH maser has been found in VLBA imaging of GIO. 624— 0.385. It would be a 12a 
and 15cr detection for the = 3^ ^ 3~ and 4+ 4~ transitions, respectively, if the average line-of-sight 
magnetic field component in the region of the J = 7/2 absorption were the same as the full three-dimensional 
magnetic field strength implied by Zeeman splitting of the 1667 MHz maser. But with only one Zeeman pair 
detected through synthesis imaging, it is impossible to conclude whether or not there exists a line-of-sight 
field reversal across the source. 

The second explanation simply states that if the magnetic field is inclined at a large angle to the line of 
sight, its projection along the line of sight may be insufficient to produce a detectable Zeeman splitting in 
Stokes V, for the case where the splitting is less than the linewidth. Assuming that the magnetic field strength 
of 6.0 mG for the single Zeeman pair detected with VLBI resolution in GIO. 624— 0.385 is a typical value for 
the source as a whole, the required inclination angle of the magnetic field to produce a line-of-sight component 
less than 1.2 mG is > 78°. If the magnetic field actually is inclined 78° to the line of sight, the expected linear 
polarization fraction of the two cr-componcnts of the Zeeman pair is sin^ 78°/ (1 + cos^ 78°) = 0.92 (Goldrcich, 
Keeley, & Kwan 1973). But no linear polarization is detected in the maser (Fish 2004), ruling out a high 
inclination at the maser site, on the east side of the H II region. Linear polarization fractions consistent with 
inclination angles as large as 76° are seen on the west side of the H II region, although most masers imply 
inclination angles much smaller than this value. Incomplete spatial coverage of OH masers across the source 
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render it difficult to estimate with certainty tlie inclination angle averaged across the source, but it is almost 
certainly less than 78°. Thus, while B^^ < B = 6.0 mG in GIO. 624— 0.385, the inclination angle is likely 
not large enough such that By < 3a = 1.2 mG. Still, an inclination angle smaller than 78° in combination 
with another factor (such as a field reversal across the source) may suffice to reduce the average line-of-sight 
magnetic field strength below our detection threshold. 

The third explanation is unlikely because a higher temperature is required to populate the ^113/2, J = 
7/2 states than for either the J ~ 3/2 or J = 5/2 states. These higher temperatures likely require that the 
distribution of OH in the J = 7/2 states be peaked closer to the central heat source than for the J = 3/2 and 
J = 5/2 states, as is noted for masers in W3(0H) (Baudry & Diamond 1998). Since the density probably 
increases with decreasing radius, the density (and therefore average field strength) of the material sampled 
by OH in the J = 7/2 states should be higher than that of the lower states. There is no evidence that the 
magnetic field strengths derived from Zeeman splitting in J = 7/2 masers (Caswell 2004) are weaker than 
those derived from masers in the J = 5/2 (Baudry et al. 1997; Caswell & Vaile 1995) or J = 3/2 transitions 
(e.g.. Fish et al. 2003). Indeed, if W3(0H) is not atypical, there is reason to believe that field strengths 
deduced from OH Zeeman splitting actually increase as measured by higher excitation states in the ^113/2 
ladder (Giisten, Fiebig, & Uchida 1994; Baudry & Diamond 1998). 

The fourth explanation differs from the third in that it suggests that masers may inherently occur at 
density enhancements in the surrounding medium. This is plausible given the conditions of formation of 

interstellar OH masers. Observed maser strengths require some combination of OH enrichment and density 
enhancement (Elitzur 1992). OH masers are believed to form in the zone between the ionization and 
shock fronts (Elitzur & de Jong 1978), where instabilities lead to inhomogeneous density enhancement 
(e.g., simulations of Garcfa-Segura & Franco 1996). The density conditions under which OH masers form 
have implications which affect the physical interpretation of their phenomenology. If masers are formed 
preferentially at density enhancements, magnetic field measurements obtained at OH maser sites should 
be higher than the average magnetic field strength in the surrounding region. It would also be strong 
evidence that observed proper motions of OH masers (as in Bloemhof, Reid, & Moran 1992) are due to 
discrete material motions, not shifting coherence paths that could be unrepresentative of the motion of the 
material.^ 

It is difficult to tell conclusively whether or not the material sampled by J = 7/2 absorption is at a 
lower density than the ^113/2 masers in most of our sources. In K3— 50, the upper limit on the magnetic field 
strength in the absorption falls within the range obtained from masers. In GIO. 624— 0.385, the inclination 
of the magnetic field to the line of sight alone is unlikely to explain the discrepancy between our upper 
limit and maser field strength values. But it is quite possible that a field direction reversal occurs across 
the source, which would reduce the effective average line-of-sight strength detectable through absorption. It 
is also not possible to rule out that there exists a field direction reversal along the line of sight. However, 
magnetic field directions deduced from OH maser Zeeman splitting in massive star-forming regions show an 
overwhelming tendency to fall into one of two categories: (1) a constant line-of-sight direction throughout 
the source, and (2) a single reversal across the source with the property that a line can be drawn dividing 
the source into a region of positive magnetic field and a region of negative magnetic field (Fish 2004). A 
reversal of the field direction along the line of sight would produce projected regions of mixed magnetic field 
direction, suggesting that they are not common in massive star-forming regions. 



The persistence of maser spot shapes also present a strong case that proper motions are due to material motions (Bloemhof, 
Moran, & Reid 1996). 
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As for W3(0H), the magnetic field measurement obtained from absorption By = 3.0 ± 0.3 mG is less 
than those obtained from masers projected atop the H II region in the ^113/2, J = 3/2 transitions (average ± 
standard error of the mean: B = 5.6± 0.7 mG, from Bloemhof, Reid, & Moran 1992), J = 5/2 transitions 
(6.6±0.5 mG, from Desmurs et al. 1998), as well as the J = 7/2 transitions (9.7±0.4 mG, from this work and 
Baudry & Diamond 1998). The average magnetic field measurements from the masers projected atop the 
H II region exceed the absorption magnetic field measurement by 3.3 a, 5.8 a, and 13.0 a for the J = 3/2, 5/2, 
and 7/2 transitions, respectively. Most masers in W3(0H) have no detectable linear polarization (Garcia- 
Barreto et al. 1988), suggesting that the magnetic field in W3(0H) is oriented close to the line of sight (i.e., 
« B) and therefore that the average magnetic field sampled by J = 7/2 absorption is less than that 
sampled by mascr emission. 

The average magnetic field strength obtained from the J = 7/2 masers is greater than those obtained 
from J = 3/2 and 5/2 masers. As previously noted, the J = 7/2 masers in W3(0H) are much more tightly 
distributed near the center of the source than in either of the other masing ^113/2 transitions. Yet the 
magnetic field strength measured from J = 7/2 absorption is significantly less than the value obtained from 
J = 7/2 masers. Taken together, these results suggest that J = 7/2 masers occur only in the highest-density 
portions of the OH cloud, even though OH likely exists in the J = 7/2 excited state throughout the source. 

Two possible scenarios could explain the prevalence of high-density material necessary for excited-state 
OH maser activity. There may be a large region of higher density coincident with the distribution of J = 7/2 

and methanol masers (Moscadelli et al. 1999) but with little small-scale clumping. Alternatively, small- 
scale (f» 10^^ cm) clumping may occur throughout W3(0H) but with increased prevalence in the region were 
J = 7/2 OH masers are observed. 

Our data do not directly distinguish between these two scenarios. Our observations of Zeeman splitting 
in excited-state OH absorption allow us to measure the magnetic field (and therefore density) averaged 
over the entire source, but not the length scale on which density fluctuations occur. Nevertheless, other 
evidence suggests that small-scale density variations are responsible for excited-state OH maser activity. 
The distribution of ground-state OH masers in W3(0H) (Reid et al. 1980) and other massive star-forming 
regions (Fish 2004) show clustering on a scale of 10"'^^ cm. This is imlikely to be caused by Kolmogorov 
turbulence, which is a scale-free process. Also, ammonia observations by Reid, Myers, & Bieging (1987) 
show that while the optical depth of NH3 is fairly constant across W3(0H), the beam-filling factor decreases 
from west to east. This suggests that the number of clumps is decreasing, not the density of any one clump. 

3.4. Line Ratios 

In local thermodynamic equilibrium, the relative strengths of the f = 4+ 3~,3"'" 3~,4+ ^ 4~, 
and 3+ — > 4~ absorption lines are 1 : 27 : 35 : 1. The ratio of the = 4+ — » 4~ to 3+ ^ 3~ main lines 
ranges from 0.41 (GlO.624-0.385, 1 km s"^) or 0.97 to 1.62 in our sources, with 1.30 being the LTE value 
for optically thin lines. The satellite lines, when detected in absorption, are always enhanced relative to the 
expected LTE value for the F = 4+ ^ 4~ line and almost always enhanced relative to the F = 3+ — > 3~ 
transition, the single exception being the F = 4+ ^ 3^ line in W3(0H). 

Matthews et al. (1986) surmise that the excitation temperatures in the hyperfine lines are not equal due 
to line overlaps in the 84 /zm lines connecting the ^n3/2, J = 5/2 and ^n3/2, J = 7/2 states of OH. Viscuso 
et al. (1985) note that collisions will preferentially excite the ^n3/2, J = 7/2 negative-parity state, but that 
the 84.42 /xm line that excites the positive-parity state may be in resonance with CO emission at 84.41 fim. 
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Matthews c!t al. point out that far infrared hne overlaps will equalize populations between the F = and 
F = 3+ states, but the F = 3~ state will depopulate relative to the F = 4~ state. Hence, the F = 4+ ^ 4~ 
and 3+ — > 4~ absorption lines will have a lower excitation temperature than the F = 3+ ^ 3~ and 4+ 3~ 
lines. 

For the two cases in which both satellite lines are detected, the F = 3+ — > 4~ line is stronger than the 
F = 4+ ^ 3^ line. Additionally, the 3+ 4~ lino is detected in two sources in which the 4+ 3~ line 
is not detected. This is consistent with a lower excitation temperature for the 3+ 4^ state than for the 
4+ — > 3~ state. However, we do not find evidence that the F = 4+ ^ 4~ lines are enhanced relative to the 
3+ — > 3~ lines, as predicted by the Matthews et al. model. 

4. Conclusions 

We have detected ^113/2, J = 7/2 OH absorption toward six massive star-forming regions. Main-line 
absorption was detected in both the F = 3+ ^ 3~ and 4+ 4~ lines toward all sources. Additionally, we 
detected at least one satellite line in absorption toward four of the six sources. We believe that this is the 
first detection of the ^113/2, J = 7/2 satellite lines in interstellar sources. 

Stokes V spectra of the main lines were produced for these six sources. In the case of W3(0H), Zceman 
splitting of the F = 3+ ^ 3^ absorption results in a magnetic field measurement of _B|[ = 3.0 ± 0.3 mG, 
consistent in magnitude and sign with the 3.1 ± 0.4 mG obtained by Glisten, Fiebig, & Uchida (1994). 
W3(0H) is the only source in which a positive detection of Zeeman splitting in the J = 7/2 absorption 
has been obtained. The component of the magnetic field along the line of sight is comparable to the full 
magnetic field strength measured in J ~ 3/2 and J = 5/2 masers, suggesting that these masers do not 
preferentially form in high density, low filling factor regions where the density significantly exceeds that of 
the surrounding, non-masing OH. However, the line-of-sight component of the magnetic field determined 
from J = 7/2 absorption is much smaller than that determined from J = 7/2 masers. 

The grand question is whether the small line-of-sight magnetic field strength measured in ^n3/2, J = 7/2 
absorption necessarily implies density enhancement at maser sites. If J = 7/2 absorption occurs only where 
the J = 7/2 masers occur, then the full magnetic field strength (as deduced from the J = 7/2 masers) in this 
region ranges from 5.6 to 11.3 mG (Baudry & Diamond 1998), consistent with the range of field strengths 
determined from J = 3/2 masers at this site at the northern limb of the H II region (Bloemhof, Reid, & 
Moran 1992). Restricting consideration to the F = 3+ ^ 3~ transition, the transition in which a positive 
result for Zeeman splitting is obtained in absorption, favors the upper end of this range of field strengths, 
since the only detected maser Zeeman pair implies a full magnetic field strength of 10.3 mG (Table 2). The 
magnetic field in this region does not appear to be significantly inclined to the line of sight (Garcfa-Barrcto 
et al. 1988), suggesting that the line-of-sight magnetic field strength is comparable to the full magnetic 
field strength. Since the average magnetic field strength derived from masers in this region is 9.1 mG, three 
times the field strength measured in absorption, it might be concluded that masers are denser than the 
surrounding, non-masing region by a factor of 9, under the reasonable assumption that B ex n^/^. 

However, it is likely that J = 7/2 absorption comes from other areas in front of the H II region as 
well. Bloemhof, Reid, & Moran (1992) identify three other regions of ^H3/2, J = 3/2 masers projected 
atop the H II region in W3(0H): a central clump with magnetic fields of 6.2 and 7.1 mG, a southern 
clump with fields of 2.3 to 6.0 mG, and a western clump with a magnetic field of 1.8 mG. Observations 
of the ^11^/2, J = 9/2 transitions, with an excitation temperature of 511 K above ground, show strong 
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absorption at the northern clump as well as weaker absorption at the central and southern clumps (Baudry 
& Menten 1995). In the "^111/2, J = 3/2 (270 K above ground) transitions, significant absorption is seen 
at and between all four ^113/2, J = 3/2 maser clumps, while in the J = 5/2 (415 K above ground) 

transitions, absorption is seen mainly along a line running through the northern, central, and southern 
clumps, with weaker absorption from the western clump (Baudry ct al. 1983). There is no published map of 
^113/2, J = 7/2 (290 K above ground) absorption in W3(0H), but the distribution of OH in other comparable 
excited states suggests that strong ^113/2, J = 7/2 absorption would be seen over the majority of the western 
half of the H II region, including the sites of all four ^113/2, J = 3/2 maser clumps. This would imply that 
the density at ^113/2, J = 3/2 maser sites is about 2 to 4 times that of the non-masing OH. On the other 
hand, 7820 MHz (2Hi/2, J = 3/2, F = 2+ -> 2") and 8190 MHz (2Hi/2, J = ^I2,F = 3" ^ 3+) emission 
is seen exclusively near the northern maser clump (Baudry et al. 1983), and 6031 MHz (^113/2, J = 5/2, 
= 2~ — > 2"*") and 6035 MHz (F = 3^ ^ 3^) maser emission is strongest near the northern maser clump 
(Moran et al. 1978; Desmurs et al. 1998). In total, this suggests that while ^n3/2, J = 7/2 absorption 
occurs over most of the western half of W3(0H), the strongest contribution to the absorption most likely 
comes from the northern clump, which has the highest average magnetic field strength. The unweighted 
average magnetic field determined from OH masers in all transitions atop the H II region is 6.9 mG. Due to 
the distribution of maser spots, this effectively gives highest weight to the northern clump and lowest weight 
to the eastern clump, consistent with our arguments based on the distribution of excited-state ^Hi/2 and 
^n3/2 OH absorption and emission. Thus, the density at maser sites is likely several (w 5) times that of the 
non-masing regions in W3(0H). More precise quantitative results would require an interferometric map of 
^H3/2, J = 7/2 absorption. 

We are able to place upper limits of several milligauss on the line-of-sight component of the magnetic 
field in two other sources. The 2.9 mG limit on K3— 50 is about 50% of the full three-dimensional magnetic 
field strengths obtained from Zeeman splitting of the J = 3/2 (Fish 2004) and J = 5/2 masers (Baudry 
et al. 1997). The 1.2 mG limit on GIO. 624— 0.385 is a factor of 5 smaller than the 6.0 mG obtained from 
a J = 3/2 maser Zeeman pair (Fish 2004). It is not clear if this occurs because the masers arc in higher 
density (and hence higher magnetic field) clumps or if other effects, such as field reversals or large angles of 
the magnetic field to the line of sight, are present in portions of the source. 

We have detected F = 3+ — > 4~ and/or 4+ — > 3~ satellite line absorption in four sources. Of these 

two, absorption in the 3"'" 4~ line is always stronger, as predicted by Matthews et al. (1986). Satellite 
line absorption appears enhanced relative to the main lines, from that which would be expected in local 
thermodynamic equilibrium. 

Maser emission was also observed toward two sources. In W3(0H) we find seven pairs of LCP and RCP 
masers implying magnetic field strengths from 6.9 to 11.3 mG. We detect line components at all velocities 
where strong maser fcatiires were previously observed, although wc do not have the spatial resolution to 
separate the multiple lines detected by Baudry & Diamond (1998) using the VLBA. In addition, we find a 
pair of previously undetected masers centered at —44.85 km s~^. In ON 1, we find the maser at 14 km s~^ 
previously detected by Baudry & Desmurs (2002). We do not see the strong maser at —0.13 km s~^ that 
they did, but we find a new weak maser centered at 0.31 km s"-'^. Maser strength variability, previously 
noted in the ^H3/2, J = 7/2 masers of W3(0H) (e.g., Baudry & Diamond 1998), appears to be operating 
in ON 1 as well. 

Support for this work was provided by the NSF through award GSSP04-000I from the NRAO. We thank 
T. Minter for help in setting up the observations, J. Braatz and G. Langston for assistance in data reduction. 
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and M. Modjaz for helpful comments in the interpretation of the data. We also thank an anonymous referee 
for helpful comments. 
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Table 2. Maser Parameters 











LCP 






RCP 






Source 


Transition 




Ta 


At; 




Ta 


Av 


B 








(km s-i) 


(K) 


(km s-i) 


(km s-i) 


(K) 


(km s-i) 


(mG) 


W3(0H) 


3+- 


3- 


-42.53 


0.15 


0.32 


-42.29 


0.15 


0.29 


10.3 




4+ - 


^ 4- 


-44.94 


1.93 


0.18 


-44.75 


1.82 


0.20 


10.5 








-43.97 


4.01 


0.38 


-43.84 


0.56 


0.22 


7.0 








-43.66 


9.77 


0.24 


-43.54 


7.02 


0.19 


6.9 














-43.51 


6.22 


0.70 


a 








-43.33 


12.31 


0.15 


-43.19 


12.36 


0.17 


8.2 








-43.08 


62.65 


0.30 


-42.93 


63.27 


0.31 


8.5 








-42.46 


64.50 


0.27 


-42.26 


63.35 


0.28 


11.3 


ON 1 


4+ - 


^ 4- 


14.09 


1.11 


0.20 


14.02 


1.14 


0.21 


-3.8 








0.38'' 


0.07 


0.26 


0.24'' 


0.06 


0.30 


-8.3'' 



'^Magnetic field strength is 8.5 mG if paired with LCP maser at —43.66 km s ^. See §3.1 for details. 

''Table velocities and magnetic field assume that these are = 4+ — »■ 4~ masers. See §3.1 for 
details. 
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Table 3. Absorption Parameters 



Source 


Transition 


fLSR 


Ta 


Aw 








/I — 1 \ 

(Jtm s ) 


(K) 


/I — 1 \ 
(km s ) 




3+ 


3" 


—47.22 


—0.177 


3.75 








—45.03 


—0.406 


2.40 




3+ 


^ 4- 


—46.03 


—0.039 


2.84 




4+ 


^ 3- 


—44.71 


—0.013 


3.40 




4+ 


4- 


—48.16 


—0.156 


2.70 








—44.99 


—0.620 


3.25 


01 fl 624— fl 385 


3+ 


^ 3- 


—2.54 


—0.330 


5.57 








1.02 


—0.114 


3.80 




3+ 


•^4- 


-2.36 


-0.030 


6.64 




4+ 


^4- 


-1.86 


-0.536 


6.13 








0.99 


-0.047 


1.89 


G28.199-0.048 


3+ 


^ 3- 


92.45 


-0.027 


6.12 




4+ 


^4- 


93.82 


-0.038 


8.64 


W49 


3+ 


^ 3- 


4.34 


-0.152 


4.11 








6.22 


-0.219 


21.69 




3+ 


^4- 


5.59 


-0.107 


23.07 




4+ 


^4- 


4.27 


-0.232 


5.07 








8.23 


-0.347 


19.87 


K3-50 


3+ 


^ 3- 


-25.07 


-0.216 


5.98 








-20.10 


-0.040 


2.90 




3+ 


^4- 


-25.50 


-0.018 


7.11 




4+ 


^ 3- 


-25.08 


-0.010 


3.19 




4+ 


-»4- 


-25.02 


-0.276 


6.34 








-20.24 


-0.039 


2.14 


ON 1 


3+ 


^ 3- 


10.88 


-0.014 


7.13 




4+ 


^4- 


10.74 


-0.016 


7.34 
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Table 4. Magnetic Field Measurements from Absorption and Maser Emission Lines 



Absorption 





3+ 3- 


4+^4- 


Masers** 




Source 


B (mG) 


B (mG) 


B (mG) 


References 


W3(0H) 


3.0 ±0.3 


< 7.0^ 


5.6^ 


9.7 


1,2,3,4 


GlO.624-0.385 


< 1.5 


< 1.2 


-6.0 




5 


G28.199-0.048 


< 60 


< 67 


6.7 




6 


W49 


< 11 


< 7.8 


-4.3 




7 


K3-50 


< 2.9 


< 3.1 


-3.7- 


-7.2 


5,8 


ON 1 


< 73 


< 97 


-2.5- 


-6.0 


1,5,8 



'^Average magnetic field strengths deduced from Zeeman splitting in one or more ^113/2 OH transitions are 
presented. Note that magnetic field strengths determined from Zeeman splitting of OH masers are full three- 
dimensional values, while those obtained from absorption are sensitive to only the line-of-sight projection of 
the full magnetic field. See §3.3 for more details. 

''Contamination from strong masers and uncertainty in the strength and location of the absorption features 
precludes obtaining a lower estimate of the Zeeman splitting. 

Note. — Upper limits are 3a values. 



References. — (1) this work; (2) Bloemhof, Reid, & Moran 1992; (3) Desmurs et al. 1998; (4) Baudry & 
Diamond 1998; (5) Fish 2004; (6) Caswell & Vaile 1995; (7) Caswell 2003; (8) Baudry et al. 1997 
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Fig. 1. — Masers in the = 4+ ^ 4 transition of W3(0H). LCP emission is plotted as a solid line, and 
RCP emission is plotted as a dotted line. Maser fit parameters are given in Table 2. 



W3(0H) F = 3 - 3 

0.2 I 1 1 1 1 , 1 , 1 1 , 1 1 1 , — 




-43 -42.5 -42 -41.5 

Velocity (km s ^) 

Fig. 2. — Masers in the F = 3+ — > 3~ transition of W3(0H). See Figure 1 caption for details. 
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Fig. 4. — Masers near km s~-^ in the = 4+ ^ 4~ transition of ON 1. The LCP and RCP data are 
plotted as histograms, and the best Gaussian fits are plotted as curves. The RCP data have been shifted by 
0.2 K for clarity. 
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¥3(0H) Stokes I F = 4 - 4 



0.5 









-50 -45 -40 \ 




F = 3"^ - 4" \ 




~< < 1 < < < < 1 < < < , 1 < 





-70 -60 -50 -40 

Velocity (km s ^) 



Fig. 5. — Stokes I spectrum for the F = 4+ — > 4~ and 3+ 4~ transtitions of W3(0H). Data points are 
Banning weighted and plotted as a histrogram. The Gaussian fit parameters listed in Table 3 are plotted 
as a curve. The masers are plotted in Figure 1 with a much larger scale to show the strong emission. See 
§3.2 for details regarding the fits to the main-line absorption features. The F = 3+ ^ 4~ line is shown in 
greater detail in Figure 6. 



W3(0H) Stokes I F = 3 - 4 

0.05 I , , , , , , , 1 , , , , , , — 




-0.1 I ' ' ^ ' ' ' ' 1 ' ' ' ' ^ ' 1 

-50 -45 -40 

Velocity (km s ^) 

Fig. 6. — Stokes I spectrum for the F = 3+ — > 4~ transtition of W3(0H). See Figure 5 caption for more 
details. 
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W3(0H) 



Stokes I 



F = 3 - 3" 



-0.2 - 



< 




-0.4 - 



-50 -45 -40 

F = 4^ - 3" 



-50 -40 

Velocity (km s 



-30 

-In 



-20 



Fig. 7. — Stokes I spectrum for tlie F = 3+ ^ 3~ and 4+ 3~ transitions of W3(0H). Sec Figure 5 caption 
for more details. The maser at —42 km s~^ is plotted in Figure 2. The F = 4+ ^ 3^ line is shown in greater 
detail in Figure 8. 



W3(0H) 



0.02 



< 



-0.02 



-0.04 




-50 



-45 



Velocity (km s 



-40 



Fig. 8. — Stokes I spectrum for the F = 4+ — > 3 transition of W3(0H). See Figure 5 caption for more 
details. 



- 20 - 



GIO. 624-0. 385 Stokes I 



F = 4 - 4" 



< 



-0.2 - 



-0.4 - 



-0.6 




-20 -10 

Velocity (km s ^) 



Fig. 9. — Stokes I spectrum for the F 
caption for more details. The F = 3+ 



4+ 4" and 3+ 4' transitions of GlO.624-0.385. See Figure 5 
* 4~ line is shown in greater detail in Figure 10. 



GIO. 624-0. 385 Stokes I 



F = 3 - 4" 



0.02 



^ -0.02 



-0.04 



-0.06 




10 



-5 
Velocity (km s ^) 



Fig. 10. — Stokes I spectrum for the F = 3^ ^ 4 transition of GIO. 624— 0.385. Sec Figure 5 caption for 
more details. The turnoff at higher velocity is due to F = 4+ ^ 4~ absorption, shown in Figure 9. 
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Fig. 11. — Stokes I spectrum for the F = 3+ ^ 3^ transition of GIO. 624— 0.385. See Figure 5 caption for 
more details. The arrow indicates where f = 4+ ^ 3^ absorption would appear if at the velocity of the 
main absorption component. 




Fig. 12. — Stokes I spectrum for the F = 4+ ^ 4 transition of G28. 199-0.048. See Figure 5 caption for 
more details. 
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Fig. 13.— Stokes I spectrum for the F = 3+ ^ 3" transition of G28. 199-0.048. See Figure 5 caption for 
more details. 



W49 Stokes I F = 4 - 4 




-40 -20 20 40 

Velocity (km s ^) 

Fig. 14. — Stokes I spectrum for the F = 4+ ^ 4~ and 3+ 4~ transitions of W49. See Figure 5 caption 
for more details. 
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W49 Stokes I F = 3 - 3 




-40 -20 20 40 

Velocity (km s ^) 



Fig. 15. — Stokes I spectrum for the _F = 3+ ^ 3~ transition of W49. See Figure 5 caption for more details. 
The arrow indicates where F = 4+ — » 3~ absorption would appear if at the velocity of the main absorption 
component. 
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F = 4 - 4" 
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-0.2 - 



Velocity (km s 



Fig. 16. — Stokes I spectrum for the F = 4+ ^ 4 and 3"'" ^ 4 transitions of K3— 50. See Figure 5 caption 
for more details. The F = 3+ ^ 4~ line is shown in greater detail in Figure 17. 
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Fig. 17. — Stokes I spectrum for the F = 3+ — > 4 transition of K3— 50. See Figure 5 caption for more 
details. The turnoff at higher velocity is due to F = 4+ ^ 4~ absorption, shown in Figure 16. 
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Fig. 18. — Stokes I spectrum for the F = S"*" ^ 3 and 4+^3 transitions of K3— 50. See Figure 5 caption 
for more details. The F = 4+ ^ 3~ line is shown in greater detail in Figure 19. 
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Fig. 19. — Stokes I spectrum for the F = 4+ 3 transition of K3— 50. See Figure 5 caption for more 
details. The turnoff at lower velocity is due to F = 3"*" ^ 3~ absorption, shown in Figure 18. 
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ON 1 



Stokes I 



F = 4 - 4 
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-10 10 

Velocity (km s ^) 



30 



Fig. 20. — Stokes I spectrum for tlie F = 4+ ^ 4 transition of ON 1. See Figure 5 caption for more details. 
The masers at and 14 km are plotted in Figures 3 and 4. 
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Fig. 21. — Stokes I spectrum for the F = 4+ — > 4 transition of ON 1. See Figure 5 caption for more details. 
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Fig. 22. — Stokes V spectrum for the F = 3+ — > 3~ transition of W3(0H). The sohd Hne shows the best fit 
for the derivative of Stokes I for the main absorption feature at —45 km s~^, and the dashed Hne shows the 

best fit for both F = ^ 3~ absorption features. The fit fines correspond to a magnetic field of +3.0 mG, 
as explained in §3.3. The large value of Stokes V near —42.4 km s~^ is due to a pair of maser features, which 
have been excluded from the fit. 
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